Renewable and Sustainable Energy Reviews 36 (2014) 123-134 


Renewable and Sustainable Energy Reviews 


journal homepage: www.elsevier.com/locate/rser 


Contents lists available at ScienceDirect 


Does energy consumption contribute to climate change? Q) oa 
Evidence from major regions of the world 


Ghulam Akhmat ®!, Khalid Zaman”, Tan Shukui ®*, Faiza Sajjad ° 


* College of Public Administration Huazhong University of Science and Technology (HUST), 1037, Luoyu Road, Wuhan, PR China 
> Department of Management Sciences, COMSATS Institute of Information Technology, Abbottabad, Pakistan 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 15 November 2013 
Received in revised form 

11 April 2014 

Accepted 13 April 2014 
Available online 9 May 2014 


Keywords: 

Energy consumption 
Atmosphere 

Topography 

Water system 

Living organisms threatened 
World 


This study is a contribution to the on-going debate over whether there is a relationship between energy 
consumption and climatic change variables, although the developed regions are among the most 
energy intensive economies in the world, little attention has been paid to the features of their energy 
consumption and climatic variations. Therefore, this study empirically investigates the two variables 
dynamic relationship in five broader regions of the world i.e., South Asia, Middle East and North Africa 
(MENA region), Sub-Saharan Africa, East Asia and Pacific and the aggregate data of the World. The major 
climatic variables include atmospheric, topographic, living organisms threatened; water system and 
growth factors are used to investigate energy—climate nexus over a period of 1975-2011. The results of 
the study indicate that there exists a long-run equilibrium relationship between energy consumption 
and climatic variables which shows climatic variations due to changes in energy consumption in 
different regions of the world. Although, the causality results are mixed among regions, we do find a 
systematic pattern. The present study find evident of unidirectional causality between the electric power 
consumption and climatic factors in the World’s selected region. Sound and effective energy consump- 
tion strategy may reduce the burden of global warming situation in the world. 

© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 


energy is produced, delivered, and consumed in the different regions 
of the World. Many factors, both natural and human, can cause 


Changes in temperature, precipitation, sea level, and the fre- changes in Earth’s energy balance, i.e., changes in the greenhouse 
quency and severity of extreme events likely affect on how much effect, which affects the amount of heat retained by Earth’s atmo- 
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sphere, variations in the sun's energy reaching Earth and changes in 
the reflectivity of Earth's atmosphere and surface [1]. The World 
Meteorological Organization (WMO) suggests 30 years as a standard 
time span for defining climate of a region. Common examples of 
climate are tropical, polar, marine and Mediterranean [2]. The future 
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climate is already set over this time period and the consequences 
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cannot be ignored [3]. Climate varies from place to place, depending 
on latitude, distance to the sea, vegetation, orography and other 
factors. Similar to weather, the climate variability also occurs at 
different time scales. The long term mean of the climate at a place is 
very important for our understanding since this determines many 
factors which are useful for human living [4]. 

The worldwide energy consumption in 2006 was close to 498 
exajoules. This is equivalent to an energy convergence of 15.8TW into 
the populated regions, where energy is consumed and dissipated into 
the atmosphere as heat. Although energy consumption is sparsely 
distributed over the vast Earth surface and is only about 0.3% of the 
total energy transport to the extratropics by atmospheric and oceanic 
circulations, this anthropogenic heating could disrupt the normal 
atmospheric circulation pattern and produce a far-reaching effect on 
surface air temperature [5]. Table 1 shows the data trend of the 
variables for the ready reference. 

The increasing threat of global warming and climate change has 
focused attention on the relationship between energy consumption 
and environmental pollutants [6]. Dolsak [7] analyzes factors affecting 
countries' commitment to mitigating global climate change within the 
scope of existing international institutions. A theoretical model of 
governments’ decision making is presented and tested for 91 coun- 
tries. The empirical analysis suggests that national commitment is 
significantly affected by the national government's incentives and the 
ability to affect the global level of GHG emissions, impacted more by 
the incentives than by the ability, and not affected by the aggregate 
levels of economic benefits. Leckebusch and Ulbrich [8] examine the 
relationship between cyclones and extreme wind events over Europe 
using global as well as regional climate model simulations. Climate 
change simulations based on the Special Report on Emission Scenarios 
(SRES) A2 and B2 are used. The results reveal that changes occur in 
particular with respect to the A2 scenario for extreme cyclone 
systems, while for B2 the changes are less pronounced. Especially 
over western parts of Central Europe, the track density of extreme 
cyclones increases for A2, accompanied by a tendency towards more 
intense systems. Swim and Becker [58] examines German versus U.S. 
residents’ (predominantly students’) efforts to engage in direct and 
indirect behaviors that lessen their personal contribution to green- 
house gases. The results show that Germans are more energy 
reduction behaviors because they were more likely to endorse bio- 
spheric environmental concerns, less likely to endorse egoistic envir- 
onmental concerns, less likely to think that personal costs of energy 
reduction behaviors were important, and more likely to think ethical 
considerations were important. Sultan [9] examines the export-GDP 
nexus and electricity-GDP nexus in Mauritius for the period of 1970- 
2009. The multivariate Granger-causality analysis indicates that elec- 
tricity and exports Granger-cause economic growth in the long-run. 
Fan et al. [10] focuses on the vulnerability of electricity demand due to 
the impact of climate change. The results show that electricity 
demand of all sectors is more sensitive to the heat pressure than cold 
pressure. Different electricity demand vulnerability to climate warm- 
ing exist in different sectors, ranked from higher to lower by primary 
industry, tertiary industry, secondary industry and residential sector. 

Al-mulali et al. [11] explored the relationship between urbani- 
zation, energy consumption, and CO, emission in the MENA 
countries during the period 1980-2009. The results show that 
there was a long-run bidirectional positive relationship between 
urbanization, energy consumption, and CO, emission. However, 
the significance of the long-run relationship between urbaniza- 
tion, energy consumption, and CO, emission varied across the 
countries based on their level of income and development. 
Borgstede et al. [12] explores public opinions regarding climate 
change and mitigation options and examines how psychological 
factors determine self-reported energy-efficient behavior. The 
results of an opinion poll conducted in 2005 and 2010 are 
compared. The number of respondents favoring new technologies 


as a way to reduce emissions was substantially lower in 2010 than 
in 2005, whereas there was an increase in the number of people 
who acknowledged that lifestyle changes are necessary to counter- 
act climate changes. An analysis of the 2010 survey revealed that 
respondents with pro-environmental attitudes towards global 
warming favor significantly increased use of renewable energy 
technologies and greater engagement in energy-efficient beha- 
viors. Ke et al. [13] examines the Energy-Saving Performance 
Contract (ESPC) of an office building by applying IPMVP Option 
D in combination with the energy analysis model established for 
the building by eQUEST simulation software to calibrate energy 
consumption simulation results using actual electricity billing 
data. The results indicate that, compared to actual energy con- 
sumption, the mean bias error (MBE) and root mean square error 
(RMSE) for uncalibrated simulation results are 24.48% and 
125,050, whereas the MBE and RMSE for calibrated simulation is 
0.37% and 34,197. When lighting power density increases or 
decreases by 50%, overall energy consumption decreases by 
30.78% or increases by 31.19%, respectively. Therefore, illumination 
density has the greatest impact on energy consumption. Streimi- 
kiene and Siksnelyte [14] examine the impact of the electricity 
market regulation on generating technologies, including renew- 
able in Lithuanian and Poland. The results reveal that the main 
driving forces behind the rationale for reform, electricity reform 
characteristics, the impact of electricity market reform on elec- 
tricity prices and electricity market reform and non-reform related 
factors that have influenced investor's choice for a specific gen- 
eration technology or a technology mix. 

Akhmat et al. [15] investigate the relationship between green- 
house gas (GHG) emissions, energy mix and carbon emissions in 
the panel of 35 developed countries. The results conclude that 
electricity production from oil, gas, and coal sources increases the 
GHG emissions and air pollution in the region, however, the 
intensity is far less that through fossil fuel. In another study of 
Akhmat et al. [16] examine the cause-effect relationship between 
environmental pollutants and energy consumption in the selected 
SAARC countries, over a period of 1975-2011. The results conclude 
that energy consumption acts as an important driver to increase 
environmental pollutants in SAARC countries. Khan et al. [17] 
examines the causal relationship between energy consumption 
and greenhouse gas emission for seven largest regions of the 
World. The results show that the energy consumption Granger 
causes greenhouse gas emission but not vice versa. Mudakkar 
et al. [18] investigate the long-run relationship between economic 
growth, energy consumption, environment and natural resources 
in Pakistan by using a data from 1975 to 2011. The results support 
the conjunction of some unidirectional causality, feedback hypoth- 
esis and neutrality hypothesis between the variables in a country. 

The above discussion confirms the strong correlation between 
energy consumption and climate change variables. In the subse- 
quent sections, an effort has been made to find the empirical 
investigation on energy consumption and five broad categories of 
climatic change variables i.e., atmosphere, topography, living 
organisms threatened, water system and economic growth in the 
context of different regions of the world. 

The study divided into following sections: after introduction 
which is presented in Section 1 above, section shows data source 
and methodological framework. Results are discussed in Section 3. 
Section 4 concludes the study. 


2. Data source and methodological framework 
The present study is based on annual time series data covering 


the time period from 1975 to 2011 for South Asia, Middle East and 
North Africa (MENA), Sub-Saharan Africa, East Asia and Pacific and 


Table 1 
Data view of the variables. 


Source: World Bank [19] and ECV [20]. 


Aggregate data of the regions 


MENA 


South Asia 


Sub-Saharan Africa 


East Asia and Pacific 


World 


Indicator Name 


Electric power consumption (kWh per capita) 

CO- emissions from transport (% of total fuel combustion) 
Arable land (% of land area) 

Urban population growth (annual %) 

Agriculture, value added (% of GDP) 

Combustible renewables and waste (% of total energy) 
GDP per capita (current US$) 

Industry, value added (% of GDP) 

GDP per unit of energy use (PPP $ per kg of oil equivalent) 
Population growth (annual %) 


Electric power consumption (kWh per capita) 

CO- emissions from transport (% of total fuel combustion) 
Arable land (% of land area) 

Urban population growth (annual %) 

Agriculture, value added (% of GDP) 

Combustible renewables and waste (% of total energy) 
GDP per capita (current US$) 

Industry, value added (% of GDP) 

GDP per unit of energy use (PPP $ per kg of oil equivalent) 
Population growth (annual %) 


Electric power consumption (kWh per capita) 

CO- emissions from transport (% of total fuel combustion) 
Arable land (% of land area) 

Urban population growth (annual %) 

Agriculture, value added (% of GDP) 

Combustible renewables and waste (% of total energy) 
GDP per capita (current US$) 

Industry, value added (% of GDP) 

GDP per unit of energy use (PPP $ per kg of oil equivalent) 
Population growth (annual %) 


Electric power consumption (kWh per capita) 

CO, emissions from transport (% of total fuel combustion) 
Arable land (% of land area) 

Urban population growth (annual %) 

Agriculture, value added (% of GDP) 

Combustible renewables and waste (% of total energy) 
GDP per capita (current US$) 

Industry, value added (% of GDP) 

GDP per unit of energy use (PPP $ per kg of oil equivalent) 
Population growth (annual %) 


Electric power consumption (kWh per capita) 

CO, emissions from transport (% of total fuel combustion) 
Arable land (% of land area) 

Urban population growth (annual %) 

Agriculture, value added (% of GDP) 

Combustible renewables and waste (% of total energy) 
GDP per capita (current US$) 

Industry, value added (% of GDP) 

GDP per unit of energy use (PPP $ per kg of oil equivalent) 
Population growth (annual %) 


1975 


385.9589 
22.31333 
4.631538 
4.418982 
6.795582 
1.912398 
1171.182 
61.61465 
4.2012 
2.818766 


102.7706 
21.24537 
42.436 
4.0256 
37.07562 
60.53887 
171.3853 
22.32644 
1.2412 
2.345947 


394.2142 
20.03103 
6.202102 
4.841381 
19.74935 
61.11137 
405.8909 
33.05886 
1.4256 
2.747709 


571.6852 
11.95853 
8.62947 

2.622882 
14.95414 
22.8652 

642.0526 
39.31555 
1.20125 

1.915889 


1354.412 
19.76923 
10.23455 
2.535521 
9.242602 
12.36979 
1447.128 
36.61383 
1.4256 
1.869887 


1980 


658.9411 

23.12628 
4.391027 
4.537251 
7.215304 
1.248463 
2536.302 
56.64745 
4.175139 

3.199124 


125.3677 
20.7176 

42.53245 
4.328923 
34.53766 
58.04291 
264.1109 
24.1388 

1.561842 
2.417817 


478.454 
19.64915 
6.324854 
4.783142 
18.19276 
58.07321 
709.7386 
38.31449 
1.5248 
2.857382 


712.8834 
12.33969 
8.782382 
2.919798 
12.2781 
20.10064 
1155.151 
40.89791 
1.546683 
1.474842 


1586.296 
19.74143 
10.35776 
2.619081 
7.569893 
11.91057 
2502.737 
36.96208 
2.009343 
1.747188 


1985 


948.1167 
26.44838 
4.631876 
4.401154 
9.013639 
1.000851 
2344.489 
42.5561 
3.909541 
3.207222 


172.4472 
17.58341 
42.63853 
3.522284 
30.8541 
50.62459 
294.9814 
25.00145 
2.097182 
2.309294 


529.9852 
18.53044 
6.50138 

4.489046 
21.34408 
55.88297 
478.9335 
34.01215 
1.640745 
2.887301 


836.9744 
12.14863 
10.02258 
3.550693 
11.1914 

19.15577 
1391.859 
38.94585 
2.294993 
1.507991 


1718.845 
20.41351 
10.6959 
2.643881 
6.946434 
12.48581 
2611.004 
34.78518 
2.768721 
1.731411 


1990 


1167.62 
24.69141 
4.777036 
4.396551 
10.17364 
0.860941 
2105.601 
43.40379 
4.100373 
3.389832 


242.837 

15.38718 
42.67473 
3.270318 
28.88884 
43.90287 
362.0093 
25.93772 
2.586965 
2.20608 


536.8865 
18.0518 
6.706573 
4.539593 
20.73675 
57.29198 
596.8458 
34.28444 
1.989211 
2.814934 


1075.381 
12.52186 
10.26451 
3.427731 
9.84212 

16.69232 
2581.006 
38.18083 
2.820867 
1.498695 


2120.54 
19.51218 
10.81553 
2.598245 
6.173187 
10.15911 
4214.156 
33.01016 
3.035617 
1.719234 


1995 


1415.052 
22.99968 
5.092618 
3.754324 
9.5155 

0.736994 
2437.975 
42.10966 
4.105525 
2.606495 


320.395 
15.16276 
42.35675 
2.918332 
26.28724 
38.25455 
385.2373 
26.76565 
3.060645 
2.022052 


511.778 
18.65647 
7.158539 
4.244026 
19.96077 
58.12695 
563.5924 
32.44271 
2.093796 
2.72972 


1368.9 
12.90185 
9.677466 
3.215836 
7.506287 
14.25618 
4194.546 
36.9917 
3.404429 
1.2042 


2195.71 
20.38437 
10.72625 
2.354013 
4.942842 
10.32191 
5303.617 
30.96061 
3.67498 
1.489475 


2000 


1746.85 
23.654 
4.760015 
2.6188 
8.151535 
0.688765 
3020.422 
45.85242 
4.462201 
1.926258 


352.2685 
11.71464 
42.61973 
2.662588 
23.55163 
34.50197 
450.7147 
25.661 
3.665799 
1.773379 


523.381 
20.22884 
7.338453 
3.905932 
17.43973 
59.22328 
516.1671 
34.08099 
2.36173 
2.698499 


1647.687 
14.29988 
10.03266 
2.947755 
5.853674 
12.81862 
3953.912 
35.08486 
3.989368 
0.95198 


2385.443 
21.27432 

10.63708 
2.181981 

3.983701 
10.06772 
5387.443 
29.39436 
4.421164 
1.324237 


2005 


2140.2 
23.31593 
4.870516 
2.764785 
6.998349 
0.61629 
4298.73 
50.65324 
4.648834 
2.020529 


424.2892 
10.61348 
41.88003 
2.648757 
19.2257 
31.36781 
696.8783 
27.92157 
4.743589 
1.536293 


543.6952 
22.11138 
8.109428 
3.822641 
17.36612 
57.65758 
861.1318 
32.65757 
2.941988 
2.6582 


2243.316 
12.2285 
10.00454 
3.0883 
4.735165 
10.10244 
4699.189 
33.89531 
4.411694 
0.756562 


2663.017 
20.41652 
10.72657 
2.236926 
3.376014 
9.589217 
7139.987 
28.427 

5.190933 
1.213446 


2011 


2704.733 
22.79868 
4.813974 
2.425913 
7.2415 
0.525761 
8462.2 
55.6235 
5.338129 
1.929511 


605.1593 
11.43566 
41.31778 

2.570275 
18.81664 
26.70122 
1416.131 

26.60109 
6.061074 
1.340424 


534.9305 
23.7628 

8.616258 
3.927188 
15.62852 
57.648 13 
1438.263 
32.27388 
3.802452 
2.708662 


3263.545 
11.60185 

9.809014 
2.601056 
4.296744 
8.240218 
8500.775 
32.86902 
5.285875 
0.689737 


3044.221 
19.42624 
10.7675 
2.11607 
3.120929 
9.813498 
10195.81 
26.71195 
6.449624 
1.179687 
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World aggregate data. The data divides in to six broad categories of 
the variables i.e., energy consumption which is measured by 
electric power consumption, atmospheric variables include aver- 
age precipitation, carbon dioxide emissions, industrial nitrous 
oxide emissions, methane emissions from livestock and drought, 
floods and extreme temperature, topographic variables includes 
urbanization, forest depletion, mineral resource depletion and 
arable land, living organisms threatened includes bird species, 
mammal species, plant species and fish species threatened, water 
system variables includes agriculture irrigated land, ground water 
discharge and combustible waste, and finally, growth variables 
includes industry value added, GDP per capita, GDP per unit use of 
energy and population growth. The data set of the variables is 
taken from World Development Indicator which is published by 
World Bank [19] and Essential Climate Variables (ECV) [20], data 
access matrix which is published by Global Climate Observing 
Systems (GCOS). 


2.1. Theoretical frame work 


The role of energy in climate change variables are underlined in 
a number of studies (see for example, Dale [21], McMichael et al. 
[22], Harry and Morad [23] and Holmes and Reinhart [24]). 
Climate is the result of a delicate balance between several 
elements, which include: 


© atmosphere, 

© water systems, 

e living organisms and 

© topography. 

These elements determine various factors that govern the 

climate. These factors are called climatic variables. Of these, the 

most important factors that are taken into account are rain, 

atmospheric pressure, wind, humidity, and temperature [25]. 
Scientists agree that humans very likely bring about changes in 

climate through various activities. This can happen at an individual 

level or at a group level. Each individual in today's world plays a 

role, directly or indirectly, in contributing his or her bit to climate 

change [26]. For example, 


è Electric power is generated mainly by thermal power plants. 
These release a huge amount of greenhouse gases (GHGs) [27]. 

è Vehicles run on petrol or diesel, both are fossil fuels. These 
release a huge amount of GHGs [28]. 

© The more people consume luxury goods, essentials, and house- 
hold goods, the more industry flourishes. Most industries are 
run on power generated from fossil fuels. To add to this, the 
more we consume the more waste we generate [29]. 

e A great deal of waste that we generate, such as plastics, does 
not get degraded and remains in the environment for many 
years, causing damage. The use of trees in large quantities for 
construction of houses leads to further depletion of forests [30]. 

e Changes in aerosol level causes scattering of sunlight and hence 
impacts global energy balance [31]. 


Climatologists examine weather patterns and other data from 
the past in order to make predictions about Earth's overall climate 
in the future. They make such predictions by creating climate 
models. Climate models are sets of mathematical equations. These 
mathematical equations describe the relationships among many 
atmospheric properties, or variables. Such properties include wind 
speed, air pressure, temperature, precipitation, and the concentra- 
tions of many gases in the air. Climate models allow scientists to 
predict how a change in each variable will affect the other 
variables and ultimately the global climate on Earth [32]. 


This study developed a comprehensive model based upon 
number of previous studies on the same topic which comprises 
six broad categories of energy consumption and climate change 
variables in the context of five broad regions of the world i.e., 
South Asia, MENA, Sub-Saharan Africa, East Asia and Pacific and 
aggregate data of World i.e. 


In (EPC), = @ +a In (ATM),+ @ In (TOPO), 
+a In (LOT),+a@4 In (WS),+ a5 In(GROWTH)+& (1) 


where ao represents intercept, @,;-a5 are the slope of the respec- 
tive variables, In(EPC) is the natural logarithm of electric power 
consumption, In(ATM) is the natural log of atmospheric variables 
i.e., average precipitation, carbon dioxide emissions, industrial 
nitrous oxide emissions, methane emissions from livestock and 
droughts, floods and extreme temperatures, In(TOPO) is the 
natural log of topographic variables i.e., forest depletion, mineral 
resource depletion, arable land, and urbanization, In(LOT) is the 
natural log of living organisms threatened i.e., bird species, 
mammal species, plant species, and fish species threatened, In 
(WS) is the natural log of water system which contains ground 
water recharge, agriculture irrigated land and combustible waste, 
In(GROWTH) is the natural log of growth variables which com- 
prises GDP per capita, industry value added, GDP per unit use of 
energy and population growth and e is the white noise error term. 
The dependent and independent variables used in this study are 
listed in Table 2. Electric power consumption (EPC) is used as a 
dependent variable for the study while, independent variables are 
atmosphere, topography, living organisms threatened, water sys- 
tem and economic growth variables. 

A simple theoretical model has been extracted from previous 
studies to show energy consumption and climate change variables 
in five broad regions of the world are shown in Fig. 1. 


2.2. Econometric framework 


Comparable to all other techniques, that utilize time series 
data, it is essential to distinguish that unless the diagnostic tools 
used account for the dynamics of the link within a sequential 
‘causal’ framework, the intricacy of the interrelationships involved 
may not be fully confined. For this rationale, there is a condition 
for utilizing the advances in time-series version. The following 
sequential procedures are adopted as part of methodology used. 

In order to confirm the degree, these series split univariate 
integration properties; we execute unit-root stationarity tests. The 
Augmented Dickey-Fuller (ADF) and Phillips Perron (PP) test is a 
suitable testing procedure that is based on the null hypothesis that 
a unit root exists in the autoregressive representation of the time 
series. 

The most common procedure in choosing the optimal lag 
length is to estimate a VAR model including all our variables in 
non-differenced data. This VAR model should be estimated for a 
large number of lags, then reducing down by re-estimating the 
model for one lag less until we reach zero lags. In each of these 
models, we inspect the values of Akike Information Criteria (AIC) 
and the Schwarz Information Criteria (SIC) criteria. The model that 
minimizes the AIC and the SBC is selected as the one with the 
optimal lag length. 

In general, five distinct models can be considered. Although the 
first and the fifth model are not that realistic and they are also 
implausible in terms of economic theory, therefore, the problem 
reduces to a choice of one of the three remaining models (models 
2-4). 


Model 1: No intercept or trend in CE or VAR. 
Model 2: Intercept (no trend) in CE, no intercept or trend 
in VAR. 


Table 2 
List of variables. 


Variables 


Dependent variable: energy consumption 


Electric power consumption (EPC) 


Atmosphere (ATM) 

Average Precipitation in depth (AP) 
CO, emissions from transport (CO2) 
Industrial Nitrous oxide emissions (INE) 
Methane emissions from livestock (ME) 


Droughts, floods and extreme temperature (DFET) 


Topography (TOPO) 

Forest depletion (FD) 

Mineral resource depletion (MRD) 
Arable land (AL) 

Urbanization (URB) 


Living organisms threatened (LOT) 
Bird species (BS) 

Mammal species (MS) 

Plant species (PS) 

Fish species (FS) 


Water system (WS) 

Ground water recharge (GWR) 
Agriculture irrigated land (AIL) 
Combustible renewables and waste (CW) 


Economic growth (GROWTH) 

GDP per capita (GDPPC) 

Industry value added (IVA) 

GDP per unit use of energy (GDPENRG) 
Population growth (PG) 


Bird Species 


Mammal Species 


Plant Species 


Fish Species | 


Agriculture Irrigated 
Land 


Combustible Waste 


Living Organisms 
Threatened 
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Measurement Expected sign Data source 
Kilo watt hour (kWh) per capita World Bank [19] 

Millimeter (mm) per year Positive ECV [20] 
% of total fuel combustion Positive World Bank [19] 
thousand metric tons of CO, equivalent Positive World Bank [19] 
thousand metric tons of CO2 equivalent Positive World Bank [19] 
% of population, average 1980-2009 Positive ECV [20] 
Adjusted savings: net forest depletion (current US$) Positive World Bank [19] 
Adjusted savings: mineral depletion (Current US$) Positive World Bank [19] 
% of land area Positive World Bank [19] 
Urban population growth (annual %) Positive World Bank [19] 
In numbers Positive ECV [20] 
In numbers Positive ECV [20] 
In numbers Positive ECV [20] 
In numbers Positive ECV [20] 
% of total freshwater withdrawal Positive ECV [20] 
% of total agricultural land Positive World Bank [19] 
% of total energy Positive World Bank [19] 
Current US$ Positive World Bank [19] 
% of GDP Positive World Bank [19] 
PPP $ per kg of oil equivalent Positive World Bank [19] 
Annual % Positive World Bank [19] 


Average 


Mineral Resource Precipitation 


Depletion 


Urbanization 


Topography 


Arable Land 


Forest Depletion Carbon dioxide 
emissions from 


transport 


Atmosphere 


Industrial Nitrous 
\ Oxide emissions 
Methane Emissions 
from Livestock 


Droughts, Floods 
and Extreme 
Temperature 


Electric Power 
Consumption 


Industry Value 
Added 


— = 


GDP per capita 


Economic Growth 
Ground Water 


Rechrage 


Population Growth 


GDP per unit use of 
\ Energy 


Fig. 1. Research framework for energy and climate change variables. 


Source: Self extract. 
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Model 3: Intercept in CE and VAR, no trends in CE and VAR. 
Model 4: Intercept in CE and VAR, linear trend in CE, no trend 
in VAR. 

Model 5: Intercept and quadratic trend in the CE intercept and 
linear trend in VAR. 


For the intention of investigating the long-run relationship 
among the variables, they must be co-integrated. In the multi- 
variate case, if the I(1) variables are linked by more than one co- 
integrating vector, the Engle-Granger procedure is not applicable. 
The test for co-integration used here is the likelihood ratio put 
forward by Johansen and Juselius [33], indicating that the max- 
imum likelihood method is more appropriate in a multivariate 
system. Therefore, this method is used in this study to identify the 
number of co-integrated vectors in the model. The Johansen and 
Juselius method has been developed in part by the literature 
available in the field and reduced rank regression. The co- 
integrating vector ‘r’ is defined by Johansen as the maximum 
Eigen-value and trace test. There is ‘r or more co-integrating 
vectors. Johansen and Juselius [33] and Johansen [34] propose that 
the multivariate co-integration methodology can be defined as 


In(EPC) = In(ATM, TOPO, LOT, WS, GROWTH) (2) 


which is a vector of P=5 elements. Considering the following 
autoregressive representation: 


K 
Ye= no+ DY TiYt-1+he 
tst 


Johansen's method involves the estimation of the above equa- 
tion by the maximum likelihood technique, and the testing of the 
hypothesis Ho, (æ = YE) of ‘r’ co-integrating relationships, where 
‘r is the rank or the matrix 2(0zrzP), ¥ is the matrix of weights 
with which the variable enters co-integrating relationships and & 
is the matrix of co-integrating vectors. The null hypothesis of non- 
cointegration among variables is rejected when the estimated 
likelihood test statistic ;{ = —nyP_,, ind —Aj} exceeds its cri- 
tical value. Given estimates of the Eigen-value (A;) the Eigen-vector 
(éi) and the weights (¥;), we can find out whether or not the 
variables in the vector are co-integrated in one or more long-run 
relationships among the dependent variables. 

If the time series are integrated at first difference, then one 
could run regressions in their first differences. However, by taking 
first differences, we drop the long-run correlation that is stored in 
the data. This means that one needs to use variables in levels as 
well. Error Correction Model (ECM) incorporates variables both in 
their levels and first differences. ECM depicts the short-run 
disequilibrium as well as the long-run equilibrium adjustments 
between variables. ECM term having negative sign and value 
between “0 and 1” specifies convergence of the model towards 
long-run equilibrium. 

Following Pesaran et al. [35], we assemble the vector autore- 
gression (VAR) of order p, denoted VAR (p), for the following 
growth function: 


p 
Ze=U+ Ù PiZt-it+Et (3) 
iZi 


where z+ is the vector of both x, and y, where y+ is the dependent 
variable defined as electric power consumption (EPC), x; is the 
vector matrix which represents a set of explanatory variables and t 
is a time or trend variable. According to Pesaran et al. [35], yy must 
be I(1) variable, but the regressor x, can be either I(0) or I(1). We 
further developed a vector error correction model (VECM) as 


follows: 
p-i p-1 

AZ = u+ at+AZzt 1+ DY YAY it È yAxe-iteét (4) 
i=1 i=l 


where A is the first-difference operator. The long-run multiplier 
matrix À as 


AyyA 
i-| YY i 


Axy Axx 


The diagonal elements of the matrix are unrestricted, so the 
selected series can be either I(0) or I(1). If Ayy = 0, then Y is I(1). In 
contrast, if Ayy < 0, then Y is I(0). 

The VECM procedures described above are imperative in the 
testing of at most one cointegrating vector between dependent 
variable y, and a set of regressors x;. To derive model, we followed 
the postulations made by Pesaran et al. [35] in Case III, that is, 
unrestricted intercepts and no trends. After imposing the restric- 
tions Ayy =0, 440 and a=O, the hypothetical function can be 
stated as the following unrestricted error correction model 
(UECM): 


A In (EPC), = Po +21 In (EPC),_1+/2 In (ATM); 
+£, In (TOPO),_,;+/34 In (LOT), _, 
+f, In(WS),_ 1+, In(GROWTH) 


p q 
+ >} f,Aln(EPC),_;+ È} PgA In (ATM),_; 
i=0 


i=l 


r § 
+ È fB,A1n(TOPO),_;+ È fy A In (LOT), i 
i=0 i=0 


t u 
+ È By,AIn(WS),_;+ DY B.A In (GROWTA),_;+ur 
i=0 i=0 


(5) 


where A is the first-difference operator and us is a white-noise 
disturbance term. Eq. (5) also can be viewed as an ARDL of order 
(p, q, T, S, t, u). Eq. (5) indicates that energy consumption tends to 
be influenced and explained by its past values. The structural lags 
are established by using minimum Akaike's information criteria 
(AIC), Pesaran et al. [35] suggested using the standard joint 
significance F test on the lagged levels variables. After regression 
of Eq. (5), the Wald test (F-statistic) was computed to differentiate 
the long-run relationship between the concerned variables. The 
Wald test can be carry out by imposing restrictions on the 
estimated long-run coefficients of EPC, ATM, TOPO, LOT, WS and 
GROWTH. The null and alternative hypotheses are as follows: 

Ho : B, = P2 = P3 = P4 = P5 = Ps = 0 (no long-run relationship) 

Against the alternative hypothesis 

Ha: P1 # P2 # P3 # P4 # Ps #Pe #0 (a long-run relationship 
exists) 

The computed F-statistic value will be evaluated with the 
critical values tabulated in Table CI (iii) of Pesaran et al. [35]. 
According to these authors, the lower bound critical values 
assumed that the explanatory variables x; are integrated of order 
zero, or I(0), while the upper bound critical values assumed that x; 
are integrated of order one, or I(1). Therefore, if the computed 
F-statistic is smaller than the lower bound value, then the null 
hypothesis is not rejected and we conclude that there is no long-run 
relationship between energy consumption and its determinants. 
Conversely, if the computed F-statistic is greater than the upper 
bound value, then variables among have a long-run level relation- 
ship. On the other hand, if the computed F-statistic falls between the 
lower and upper bound values, then the results are inconclusive. 

The Granger causality test is then used to determine the 
direction of causality among the six broad categories of 
the variables in this study. However, if the variables are I(1) and 
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cointegrated, the Granger causality test within the first difference 
VAR model will be misleading [36]. In such circumstances, Granger 
causality should be tested using the ECM as follows: 


A In (EPC), = a, + > k1iA In (EPC), _;+ x v1iA In (ATM),_; 


i=1 


+ È 61iA In (TOPO),_;+ 2 01iA In (LOT); 


k 
+2 yliA In (WS),_ ; 


O1iA In(GROWTH) + 91e,_1+é1t (6) 


+ 
l hae N 


k k 
A In (ATM), = a, + R2 k1iA ln (ATM),_;+ 2, V1iA In (EPC), _; 


+ Š 61iA In (TOPO), _;+ > 01iA In (LOT), _; 


i=0 


+ È yliA In (WS),_;+ > 011A In(GROWTH) 


i=0 


+91e,_14+é1t (7) 


A In (TOPO), = a; + È k1iA In (TOPO), _;+ 2 V1iA In (ATM),_; 


i=1 


+ 3 61iA In (EPC), _;+ È O1iA In (LOT); 


i=0 
4 > y1iA In (WS),_;+ $ 01iA In(GROWTH) 
i=0 i=0 
+1e,_,+élt (8) 
A In (LOT), =a, + Š k1iA In (LOT) _;+ > V1iA In (ATM),_; 
t= 


+ > 61iA In (TOPO), _;+ 2 O1iA In (EPC); _; 


i=0 
+ x y1iA In (WS),_;+ ¥ 01iA In(GROWTH) 
eee (9) 
A In (WS), = a1 + > k1iA In (WS), _;+ > v1iA In (ATM), _; 
= 
+ 2 61iA In (TOPO), _;+ 2 O1iA In (LOT), _; 
+ > y1iA In (EPC),_;+ > 011A In(GROWTH) 


i=0 
+ dei +élt (10) 


A In (GROWTH), = a7 + > k1iA In (GROWTH),_ ; 


i=1 


+ $ v1iA In(ATM),_;+ > 61iA In (TOPO), _; 


i=0 
k k 

+ ¥ @1iA In(LOT),_;+ È y1iA In(Ws),_; 
i=0 i=0 


k 
+ ¥ 011A In(EPC) 
i=0 


+01e,_1+E1t (11) 


The equations above consist of short and long-run elements. A 
and In are the notations for first difference and natural logarithm, 
respectively. The residuals €1t are assumed to be normally 
distributed and white noise. From the above equations, €;_; is 
the one period lagged error-correction term derived from the 


Table 3 
The unit root test results. 


Variables ADF ADF 


ADF (Sub- ADF (East Asia ADF 
(South (MENA) Saharan Africa) and Pacific) (World) 
Asia) 
EPC (1 (1) I(1) (1) (1) 
AP (0 (1) (1) 1(0) (1) 
CO, (1 1(1) I(1) 1(0) 1(0) 
INE (1 1(0) 1(0) (1) I(1) 
ME (1 1(0) 1(0) (1) I(1) 
DFET (1 1(1) (1) 1(0) 1(0) 
FD (0 1(1) I(1) 1(0) 1(0) 
MRD (0 (1) (1) 1(0) (1) 
URB 1(0 I(1) I(1) 1(0) I(1) 
AL (0 I(1) I(1) 1(0) I(1) 
BS* 1(0 1(1) I(1) 1(0) 1(0) 
MS? (1 (1) I(1) KJ (1) 
PSs* (1 (1) (1) 1(0) (1) 
FS* (1 (1) (1) 1(0) (1) 
AIL (1 (1) (1) 1(0) (1) 
GWR (1 (1) Ki) 1(0) I(1) 
cw 1(0 I(1) 1(0) 1(0) I(1) 
GDPPC 1(0 1(0) 1(0) 1(0) 1(0) 
IAV (1 (1) (1) (1) 1(0) 
GDPENRG I(1 1(0) I(1) (1) 1(0) 
PG (1 1(0) (1) (1) 1(0) 
Note: The null hypothesis is that the series is non-stationary, or contains a unit root. 


The rejection of the null hypothesis is based on MacKinnon critical values. The lag 
length are selected based on SIC criteria, this ranges from lag zero to lag four. The 
model with constant and trend was used to estimate the unit root tests. I(1) shows 
the variables are non-stationary at level but stationary becomes after first 
difference. I(0) shows the variables are stationary at their level. 


* Where data covers the period 199-2011. 


cointegrating equation. However, in the absence of cointegration, 
the €;_; term will be excluded from the above equation, thus it 
remains as the standard first difference VAR model. To examine 
the short run Granger causality, we apply the standard Wald test 
on the first difference lagged explanatory variables. However, to 
examine the long-run Granger causality, we employ the standard 
Wald test on both the first difference lagged explanatory variables 
and the €,_; term. 


3. Results 


The standard Augmented Dickey-Fuller (ADF) and Phillips Perron 
(PP) unit root test was exercised to check the order of integration of 
these variables. The results obtained are reported in Table 3. Based on 
the ADF unit root test statistic, it was concluded that most of the 
variables are non-stationary at their level but stationary after taking 
first difference i.e., I(1) variables, however, remaining of the variables 
are stationary at their level i.e., I(0) variables. 

The relationship between dependent variable (i.e. energy 
consumption) and the independent variables (i.e., atmosphere, 
topography, living organisms threatened, water system and eco- 
nomic growth) is observed using the multivariate cointegration 
methodology proposed by Johansen [37] and Johansen and Juse- 
lius [33]. Johansen's Cointegration Test designates at least one co- 
integrating vector. Thus, long-run relationship is maintained by 
the data generating method. Using Johansen and Juselius [33] 
multivariate cointegration tests, the study finds that a statistically 
significant relationship exists between climate change variables on 
energy consumption in different regions of the world. The follow- 
ing cointegrating vector has been determined in Table 4. 

As evident in Table 4, there are three cointegration equations in 
South Asia, four in MENA region, seven cointegration equations in 
Sub-Saharan Africa, six in East Asia and Pacific region and five 
cointegration equations in the overall data of the world, involving 
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Table 4 
Cointegration test results. 


4A: Multivariate Johansen-Juselius cointegration test 

Regions South Asia MENA 

Number of cointegrating 
equations 


Three cointegration 
equations 


Model 


Four cointegration 
equations 


Regions 


4B: Bounds testing approach to cointegration 

South Asia Fepc(EPCIATM, TOPO, LOT, WS, GROWTH) 
MENA Fepc(EPCIATM, TOPO, LOT, WS, GROWTH) 
Sub-Saharan Africa Fepc(EPCIATM, TOPO, LOT, WS, GROWTH) 
East Asia and Pacific Fepc(EPCIATM, TOPO, LOT, WS, GROWTH) 
World Fepc(EPCIATM, TOPO, LOT, WS, GROWTH) 


Significance level (%) Lower bounds, I(0) 


Narayan [39] critical value bounds of the F-statistic: intercept and no trend 
1 4.590 


5 3.276 
10 2.696 

South Asia MENA 
Diagnostic tests 
Adjusted R? 0.452 0.725 
F-statistics 112.45* 78.352* 
JB Corman) p > 0.050 p > 0.050 
LM (rra) p > 0.050 p < 0.050 
ARCH (rcn) p > 0.050 p < 0.050 
Ramsey RESET (y2spr) p > 0.050 p > 0.050 


CUSUM 
CUSUM of squares 


Stable at 5% level 
Stable at 5% level 


Stable at 5% level 
Stable at 5% level 


Sub-Saharan Africa East Asia and Pacific World 
Seven cointegration 


equations 


Six cointegration 
equations 


Five cointegration 
equations 


Wald F-statistics 


11.28* 
6.96* 
13.28* 
212.85* 
172.95* 
Upper bounds, I(1) 
6.368 
4.630 
3.895 
Sub-Saharan Africa East Asia and Pacific World 
0.614 0.385 0.458 
16.852* 3.258** 4.752* 
p > 0.050 p < 0.050 p > 0.050 
p > 0.050 p > 0.050 p > 0.050 
p > 0.050 p > 0.050 p > 0.050 
p > 0.050 p < 0.050 p < 0.050 


Stable at 5% level 
Stable at 5% level 


Stable at 5% level 
Stable at 5% level 


Stable at 5% level 
Stable at 5% level 


* Denote significance at the 1% levels. 
** Denote significance at the 5% levels. 


variables i.e., EPC, ATM, TOPO, LOT, WS and GROWTH respectively. 
The presence of the cointegration vectors shows that there exists a 
long-run relationship among the variables. This study further 
examined the existence of long-run relationship between the 
variables. We used a Hendry's general-to-specific modeling 
approach and selected the maximum lag order of three for the 
conditional ARDL-VECM. Maximum lag order of three years is 
sufficient to capture the system's dynamics' for the yearly data 
analysis [38]. Table 4B reports the results of bounds test with 
F-statistics when each variable is considered as a dependent 
variable in ARDL-OLS regressions. Based on the Narayan [39] 
critical values, all model specification i.e., Fgpc is significant at 1% 
level. Thus the null hypothesis of no cointegration is rejected in all 
regions of the world, implying long-run cointegration relation- 
ships between the variables. 

A number of diagnostic tests are conducted to ensure that the 
selected autoregressive distributed lag (ARDL) model is appropri- 
ate. Overall, in majority of the region, the selected ARDL model is 
free from heteroskedasticity, serial correlation and general speci- 
fication error, except, normality issue has found in East Asia and 
Pacific region, autocorrelation problem exists in MENA region, 
heteroskedastisity problem in MENA region and specification error 
of the model in East Asia and Pacific and World. Moreover, both 
CUSUM and CUSUM of square tests reject the presence of struc- 
tural break. Table 4, panel B, presents the calculated F-statistic for 
cointegration, the critical values provided by Narayan [39] and the 
diagnostic tests. We find that the calculated F-statistic is greater 
than the 1% upper bounds critical values. Hence, we reject the null 
hypothesis of no cointegration relationship among the variables. 
This is corroborated by the results of the Johansen-Juselius 
cointegration test. Since the suggestions of both cointegration 
tests are consistent, we confirm that the variables are cointegrated 
and that the cointegration results are robust. 


In order to check the stability of the long-run and short-run 
relationship among energy and climate change variables, we 
assess the Error Correction Model in Table 5. 

The long-run results are presented in Table 5A. In general, 
majority of cases indicate that climatic variables positively affect 
electric power consumption in different regions of the world. 
However, the intensity of the coefficient varies region to region. In 
some regions, climatic variables are more elastic in nature while in 
other regions, there is less elastic in nature which implies that 
little change in electric power consumption leads to greater 
change in climatic variables or vice versa. The overall results 
conclude that climate change leads to species range shifts and 
consequently to changes in diversity [40]. Climate change affects 
ocean conditions, including temperature, salinity, ice coverage, 
currents, oxygen level, acidity and consequently growth, body size, 
distribution, productivity and abundance of marine species, 
including those that are exploited by fisheries. Over a range of 
greenhouse gas emission scenarios, changes in the marine envir- 
onment are predicted to be more rapid in the 21st century with 
implications for marine ecosystems and dependent industries [41]. 
The increased energy efficiency, implementation of energy savings 
projects, energy conservation, and energy infrastructure outsour- 
cing reduce the level of pollution produced by urban areas [11]. 

As the variables are cointegrated, the short run elasticities are 
evaluated using the ECM. Table 5B shows the result of short-run 
elasticities and diagnostic tests. The coefficient for error correction 
term is negative and statistically significant at the 1% and 5% level, 
except MENA region, where ECM term is insignificant. The result 
shows the presence of a long-run relationship among energy and 
climate variables. The result shows that range between 10.2% mini- 
mum to 56.8% maximum variation of energy consumption in South 
Asia, Sub-Saharan Africa, East Asia and Pacific and the World are due 
to disequilibrium. Moreover, if the system exposes to shock, energy 
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Table 5 
Long-run and short-run elasticities dependent variable: A In(EPC). 


Variables Coefficients (South Asia) Coefficients (MENA) 


5A: Long-run results 


AP Positive, < 1* Positive, < 1* 
CO2 Positive, < 1* Positive, > 1 
INE Negative, <1 Positive, < 1* 
ME Positive, < 1* Positive, < 1* 
DFET Positive, < 1* Negative, <1 
FD Negative, <1 Positive, < 1* 
MRD Positive, < 1* Positive, < 1* 
URB Positive, < 1* Positive, < 1* 
AL Positive, < 1* Positive, < 1* 
BS? Negative, <1 Negative, < 1 
Ms? Negative, <1 Negative, <1 
PS? Positive, <1 Positive, <1 
FS? Positive, < 1 Positive, < 1 
AIL Positive, < 1* Positive, < 1*** 
GWR Positive, < 1* Positive, <1 
Cw Positive, > 1* Negative, <1 
GDPPC Positive, < 1* Negative, <1 
IAV Positive, < 1* Positive, < 1* 
GDPENRG Positive, < 1* Negative, < 1 
PG Positive, > 1* Negative, <1 


5B: Short-run results 


Coefficients (Sub-Saharan Africa) 


Positive, < 1** 
Positive, < 1* 
Positive, < 1* 
Positive, < 1* 
Positive, > 1*** 
Positive, < 1* 
Positive, < 1* 
Positive, < 1* 
Positive, < 1* 
Positive, < 1* 
Positive, <1 
Positive, <1 
Negative, <1 
Negative, <1 
Positive, < 1* 
Negative, <1 
Positive, < 1* 
Positive, > 1* 
Positive, < 1* 
Positive, < 1* 


Coefficients 
(East Asia and Pacific) 


Coefficients (World) 


Positive, <1 Positive, <1 

Positive, > 1** Positive, < 1* 
Positive, < 1* Positive, > 1* 
Positive, < 1* Positive, < 1* 
Positive, < 1* Positive, < 1* 
Negative, <1 Positive, < 1* 
Positive, > 1** Positive, < 1* 
Positive, < 1* Negative, <1 
Positive, < 1* Negative, <1 
Negative, <1 Positive, < 1* 
Positive, < 1* Positive, < 1* 
Positive, < 1* Positive, < 1* 
Positive, < 1* Positive, <1 

Positive, < 1* Positive, <1 

Positive, < 1* Positive, <1 

Positive, < 1* Positive, < 1* 
Positive, < 1* Positive, < 1* 
Positive, < 1* Positive, <1 

Positive, < 1* Positive, < 1* 
Positive, <1 Positive, <1 


AP Positive, <1 Positive, <1* Positive, < 1* Positive, <1 Positive, < 1* 
COz Positive, < 1* Positive, <1 Positive, < 1* Negative, > 1* Positive, <1 
INE Positive, < 1* Negative, <1 Positive, < 1* Positive, < 1* Positive, < 1* 
ME Positive, > 1* Positive, <1 Positive, < 1* Positive, < 1* Positive, < 1* 
DFET Positive, < 1* Positive, < 1* Positive, <1 Negative, < 1* Positive, <1 
FD Positive, < 1* Negative, < 1* Positive, < 1* Positive, > 1* Positive, < 1* 
MRD Positive, < 1* Negative, < 1* Positive, < 1* Positive, <1 Positive, < 1* 
URB Positive, < 1* Positive, <1 Negative, < 1* Positive, <1 Positive, < 1* 
AL Positive, < 1* Positive, <1 Positive, > 1* Positive, < 1* Negative, < 1* 
BS? Negative, <1 Positive, < 1* Positive, < 1 Positive, < 1 Positive, < 1* 
MS? Negative, <1 Positive, < 1* Positive, <1 Positive, <1 Positive, < 1* 
PS? Positive, <1 Positive, < 1* Negative, <1 Positive, <1 Positive, <1 
FS? Positive, <1 Negative, <1 Positive, <1 Positive, <1 Positive, < 1* 
AIL Positive, < 1* Positive, <1* Positive, <1 Positive, < 1* Positive, <1 
GWR Positive, > 1* Positive, <1 Negative, < 1* Positive, < 1* Negative, < 1* 
CW Positive, < 1* Negative, <1 Positive, < 1* Positive, <1 Positive, < 1* 
GDPPC Negative, < 1* Positive, <1 Positive, < 1* Positive, <1 Positive, < 1* 
IAV Positive, < 1* Positive, < 1 Negative, > 1* Positive, < 1* Positive, < 1* 
GDPENRG Positive, < 1* Positive, < 1* Positive, < 1* Positive, > 1* Positive, > 1* 
PG Positive, > 1* Negative, <1 Positive, > 1* Positive, <1 Positive, < 1* 
€t-1 —0.485** 0.314 —0.102** —0.568* —0.274* 
Diagnostic Tests 

Adjusted R? 0.401 0.565 0.587 0.412 0.389 
F-statistics 68.25* 118.58* 6.958* 2.258** 5.258* 

JB horman) p > 0.050 p > 0.050 p > 0.050 p < 0.050 p > 0.050 

LM (erran) p > 0.050 p < 0.050 p < 0.050 p > 0.050 p < 0.050 
ARCH (73 acy) p > 0.050 p < 0.050 p > 0.050 p > 0.050 p > 0.050 
Ramsey RESET (pcr) P > 0.050 p <0.050 p > 0.050 p < 0.050 p < 0.050 


CUSUM 
CUSUM of squares 


Stable at 5% level 
Stable at 5% level 


Stable at 5% level 
Stable at 5% level 


Stable at 5% level 
Stable at 5% level 


Stable at 5% level 
Stable at 5% level 


Stable at 5% level 
Stable at 5% level 


Probability values are quoted in square brackets. LM (I) tests for the null of 1st order serial correlation amongst the residuals; Het: a test based on regression of squared 
residuals on a constant and squares of the fitted values; ARCH: a test for first-order autoregressive conditional Heteroscedasticity effects; RESET: Ramsey's Regression 
Specification Error/'-test with (m, n) degrees of freedom; and the Jarque-Bera X2(2) LM test for normality of residuals. 


man 
values of the data. 


* Denote statistical significance at the 1% levels. 
** Denote statistical significance at the 5% level. 
** Denote statistical significance at the 10% levels. 


consumption takes nearly a year to converge to its long-run equili- 
brium. Regarding all other explanatory variables, we find that in the 
short run, climatic variables by and large are statistical significant 
at the 1%, 5% and 10% level. The value of R? adjusted indicates the 
ranges between 38.9% and 58.7% variation in dependent variable has 
been explained by variations in independent variables. F-value is 
higher than its critical value suggesting a good overall significance of 
the estimated model. Therefore, fitness of the model is acceptable 


indicates that the continous time serial data is not available for the said variables, therefore, forward and backward interpolation technique is used to fill the missing 


empirically. The model seems to be robust to various departures from 
standard regression assumptions in terms of residual correlation, like 
Heteroscedasticity, Autoregressive Conditional Heteroscedasticity 
(ARCH), misspecification of functional form, or non-normality of 
residuals. This result tends to suggest that, in majority of the world 
regions, the impact of any structural change over the entire sample 
period does not appear to be significant at least in terms of model 
stability. 


132 G. Akhmat et al. / Renewable and Sustainable Energy Reviews 36 (2014) 123-134 


Table 6 
Long-run Granger causality test - VECM (Causality between energy consumption 
and climate change variables). 


Variables Causality Causality Causality Causality Causality 
(South (MENA) (Sub- (East Asia (World) 
Asia) Saharan and Pacific) 
Africa) 
Energy consumption to atmosphere 
EPC to AP Jv - v v - 
EPC to CO2 v v v v v 
EPC to INE v v v - - 
EPC to ME v - v v v 
EPC to DFET v v - v v 
Energy consumption to topography 
EPC to FD v v v - v 
EPC to MRD v - - - - 
EPC to AL - v v v - 
EPC to URB v - v - v 
Energy consumption to living organisms threatened 
EPC to BS - - - - - 
EPC to MS - - - - - 
EPC to PS - v - - v 
EPC to FS v - - v - 
Energy consumption to water system 
EPC to GWR v v - - v 
EPC to AIL v - v v v 
EPC to CW v v - - - 
Energy consumption to economic growth 
EPC to GDPPC v v - v v 
EPC to IVA - v - v v 
EPC to v v v - - 
GDPENRG 
EPC to PG v - v - v 
Atmosphere to energy consumption 
AP to EPC - - - - - 
COz to EPC - - - - - 
INE to EPC v v v v v 
ME to EPC - - - - - 
DFET to EPC v - - v v 
Topography to energy consumption 
FD to EPC - - - - - 
MRD to EPC - - - - - 
AL to EPC - - - - - 
URB to EPC v v v v v 
Living organisms threatened to energy consumption 
BS to EPC - - - - - 
MS to EPC - - - - - 
PS to EPC - - - - - 
FS to EPC - - - - - 
Water system to energy consumption 
GWR to EPC) - - - - - 
AIL to EPC - - - - - 
CE to EPC - - - - - 
Economic growth to energy consumption 
GDPPC to EPC v v - v v 
IVA to EPC v v v - v 
GDPENRG to v v - v - 
EPC 
PG to EPC v v v v v 


Note: v indicates causality exists between the variables while dash (-) represents 
no causality between the variables. 


As the variables are cointegrated, we proceed to examine the 
short and the long-run Granger causality in the ECM frame- 
work. Table 6 presents the results of Granger causality among 
energy consumption and climate change variables. The overall 
results conclude that at large, there is unidirectional causality 
running towards electric power consumption to climate change 
variables which implies that electric power consumption 


changes the global weather, sound and effective energy con- 
sumption strategy may reduce the burden of global warming 
situation in the world. 

The long-run causalities have some important policy implica- 
tions in the context of different regions of the world. First, the 
finding of unidirectional Granger causality between energy con- 
sumption and climate change variables suggests that per capita 
electric power consumption is the stimulating input for enhancing 
climatic variables and economic growth [42]. The findings further 
indicate the neutrality hypothesis which treats that energy consump- 
tion as an insignicant part of climate change variables and economic 
output which shows no causality between these variables [43]. 

The awareness of global environmental challenges makes it 
essential that there is some understanding of the causal effects of 
energy consumption on climate change and development. The 
result shows that there is a bidirectional causality running 
between industrial nitrous oxide emissions to electric power 
consumption in South Asia, MENA and Sub-Saharan African region. 
Similarly, there is a bidirectional causality running between 
droughts, floods and extreme temperature to electric power 
consumption in South Asia, East Asia and Pacific and the World. 
Among topographic variables, only urbanization is the only vari- 
able which cause electric power consumption and vice versa in 
South Asia, Sub-Saharan Africa and World. Most of the growth 
variables have shown bidirectional causality with the electric 
power consumption in different regions of the world which 
implies that any policy to reduce energy consumption (and hence 
emissions) would have an effect on climate change variables and 
overall region's development. 


4. Conclusion 


This study presents new theoretical and empirical insights into 
the analysis of the causal relationship between energy consump- 
tion and climate change variables when considering a sample of 
five broader regions of the world. The results support four 
alternative and equal plausible hypotheses among different vari- 
ables in the different regions of the world. First, the “conservation 
hypothesis” that shows a unidirectional causality running from 
climatic variables to energy consumption. The energy conservation 
policies resulting in a reduction of energy consumption have little 
or no negative effects on climatic variables i.e., industrial nitrous 
oxide emission in East Asia and Pacific region and World, urbani- 
zation in MENA and East Asia and Pacific region, industry value 
added in South Asia and Sub-Saharan Africa, GDP per energy use 
in East Asia and Pacific region, and finally, population growth in 
East Asia and Pacific region. Secondly, the “growth hypothesis” 
which shows climatic change is dependent on energy use, which 
means there is a unidirectional causality running from energy to 
climatic factors. A significant change in climatic factors, knowing 
that energy use is not positive, it is due to improper use or to 
energy inefficiency i.e. which leads to distortions in climatic 
factors. Thirdly, the “feedback hypothesis” that suggests the 
existence of bidirectional causality between energy consumption 
and climatic factors in different regions of the world. In this case 
an energy policy oriented towards improving the energy efficiency 
would little or no harm climatic factors. Finally, the “neutrality 
hypothesis” that shows there is no causal relationship. Similar to 
the assumption of conservation hypothesis, an energy conserva- 
tion policy has no or little effect on climatic factors. 

The following conclusion has been emerged from the analysis i.e. 


è The climate change poses challenges to energy suppliers and 
users in terms of potential increase in the economic cost of 
energy use [44]. 
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© Rise in the average and extreme temperatures increases elec- 
tricity demand for cooling down temperatures in different 
regions of the world [45]. 

e Climate change has larger impact on the peak electricity 
demand than on the average monthly electricity demand [46]. 

© Climate change risks oil and gas supply activities in vulnerable 
coastal areas, offshore production areas, and tundra areas [47]. 

© Both climate change and rising concentrations of atmospheric 
carbon dioxide affect bioenergy production potentials [48]. 

© Expected seasonal and chronic water scarcity present risk of 
disruptions in the supply of electricity in many different 
regions of the world [49]. 

© Climate change affects the geographical pattern of renewable 
energy supply in the different regions of the world [50]. 

è Expected reductions in precipitation in the form of snowfall 
reduce hydropower production [51]. 

© In most of the cases analyzed, adaptation measures can reduce 
risks and prospects of negative consequences of energy supply 
and use [52]. 

è Energy system resilience against climate change can be 
achieved by investing in the research and development pro- 
grams [53]. 

© Most of the vulnerabilities and risks from energy supply and 
demand reflect relatively fine-grained place-based differences 
in situations [54]. 

è The variability of risks from weather-related events in both 
time and space increases with climate change [55]. 

© Impact of Climate change interacts with regulatory environ- 
ment [56]. 

© Despite uncertainties about the impact of climate change in the 
future, robust risk management strategies can be developed in 
an iterative manner that incorporates continuous observation, 
evaluation, and learning [57]. 

è A critical step in developing RRMS (robust risk management 
strategy) is to conduct vulnerability assessment of energy users 
and suppliers in advance [58]. 

© Improvement in the existing knowledge about vulnerabilities 
and risk management strategies is essential for effective cli- 
mate change risk management in the energy sector [59]. 

è A self-sustaining long term assessment plan with commitment 
is needed to improve the scientific base for streamlining long 
term solutions of climate change [60]. 

© Self-sustaining assessment plans provide value to all stake- 
holders of climate change in the regions [61]. 


These conclusion lead us to join our heads to have the efficient 
and sustainable energy consumption which have less or almost 
negligible impact on our nature format. 
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